Holographic baryonic matter in a background magnetic field 



Florian Preis/'B Anton Rebhan/'Q and Andreas Schmitt 1 ^ 
1 Institut fur Theoretische Physik, Technische Umversitat Wien, 1040 Vienna, Austria 

We discuss the effect of baryonic matter on the zero-temperature chiral phase transition at finite 
chemical potential in the presence of a background magnetic field. The main part of our study is 
done in the deconfined geometry of the Sakai-Sugimoto model, i.e. at large N c and strong coupling, 
with non-antipodal separation of the flavor branes. We find that for not too large magnetic fields 
baryonic matter completely removes the chiral phase transition: chirally broken matter persists 
up to arbitrarily large chemical potential. At sufficiently large magnetic fields, baryonic matter 
becomes disfavored and mesonic matter is directly superseded by quark matter. In order to discuss 
the possible relevance of our results to QCD, we compute the baryon onset in a relativistic mean-field 
1 model including the anomalous magnetic moment and point out the differences to our holographic 

I , calculation. 
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I. INTRODUCTION 
A. Context 



In the hadronic phase of Quantum Chromodynamics (QCD) the approximate SU(Nf)i, x SU(Nf)n chiral symmetry 
is spontaneously broken by a chiral condensate. This condensate of quark-antiquark and quark-hole pairs is expected 
_ to melt for large temperatures T and/or large baryon chemical potentials [Ib- While the finite-temperature behavior 
■ at (is = can be understood with the help of lattice calculations [l| , there are currently no first-principle methods to 
understand the fate of the chiral condensate at small temperatures and large (1b- We do know from first principles that 
. chiral symmetry is also spontaneously broken at asymptotically large [Ib - albeit by a different mechanism, a diquark 



condensate in the color-flavor locked (CFL) phase Studies of dense matter at large, but not asymptotically 

, large, (Ib currently rely on model calculations. The logical possibilities regarding chiral symmetry, borne out in these 
<^> models, are either an intermediate chirally restored region, i.e., at least two phase transitions between ordinary nuclear 
matter and CFL matter [1, @ or a continuous transition such that chiral symmetry is broken throughout the T = 
axis of the QCD phase diagram In this paper we shall ignore diquark condensation for simplicity and only 

>> \ ' consider chiral symmetry breaking via the "usual" mechanism. 

Phenomenologically this region of intermediate [Ib and small temperatures is relevant for the physics of neutron 
stars. The question of the chiral phase transition(s) in the QCD phase diagram is thus related to the question 
whether neutron stars are hybrid stars, i.e., whether they have a quark matter core surrounded by a mantle of nuclear 
matter, or whether they are entirely made out of nuclear (and possibly hyperonic) matter. Since neutron stars have 
large magnetic fields, it is interesting to investigate the chiral (and deconfinement) transitions in a large magnetic 
background. Many studies about the phases of QCD in a background magnetic field have recently been pursued, 
for instance on the lattice @ as well as in model calculations (lfj| - fl5j . From a simple estimate we can infer that 
extraordinarily large magnetic fields are required to have a sizable effect on QCD properties: if we simply convert 
the QCD scale into a scale for the magnetic field, we obtain Aq CD ~ (200 MeV) 2 ~ 2 x 10 18 G. This is large, but 
perhaps not out of reach for the cores of neutron stars whose surface magnetic fields have been measured to be as 
large as ~ 10 15 G possibly leading to magnetic fields in the core up to ~ 10 18 G or, in the case of quark 
stars, even up to ~ I0 20 G (f^- Besides this phenomenological motivation it is also of theoretical interest to consider 
large magnetic fields. For instance, it is instructive to consider the limit of asymptotically large magnetic fields which 
often becomes simple and which in turn may be used to understand fundamental properties of dense matter, also in 
a smaller magnetic background. 



'Electronic address: fprcis@hcp.itp.tuwien.ac.at 
^Electronic address: rcbhana@hep.itp.tuwien.ac.at 
JElectronic address: aschmitt @hep.itp.tuwien.ac. at | 



2 



B. Purpose 

In this paper, we use the Sakai-Sugimoto mode l fl9l.[20l to study the chiral phase transition in a magnetic field. This 
model, making use of the gauge/gravity duality [2l| - |24l |. is a top-down approach to large- N c QCD at low energies. In 
the original version of the model, Nf right-handed and Nf left-handed flavor branes are asymptotically separated on 
antipodal points of the circle of a compactified extra dimension. In this version the phase structure is very simple, at 
least in the probe brane approximation Nf -C N c . The non-antipodal, or "decompactified" , version, which corresponds 
to a (non-local) Nambu-Jona-Lasinio (NJL) model [25|,[26j, is more susceptible to a chemical potential and a magnetic 
field, even in the probe brane approximation. Therefore, the phase structure is very rich and may even be more 
appropriate to learn something about real-world QCD with N c = 3. 

The present work is an extension of our previous study [27j |. The main observation of ref. (27j was that, for certain 
intermediate values of the baryon chemical potential, chiral symmetry gets restored upon increasing the magnetic 
field. We have termed this effect "Inverse Magnetic Catalysis" (IMC) since it is in apparent contrast to "Magnetic 
Catalysis" (MC) [HI, [29| which, simply put, says that a magnetic field works in favor of a chiral condensate. The 
original works on MC have pointed out that at weak coupling fermion-antifermion pairing in the presence of a magnetic 
field is not unlike Cooper pairing in a superconductor [3(j, leading to a chiral condensate even for arbitrarily small 
attractive interaction. IMC at finite chemical potential can occur as follows. A nonzero chemical potential creates an 
asymmetry between fermions and antifcrmions and thus imposes a stress on the chiral condensate. Eventually, chiral 
symmetry gets restored at some critical chemical potential. The simple but important observation of ref. (27j was that 
the free energy cost to maintain the chiral condensate in spite of a nonzero chemical potential is proportional to the 
magnetic field. At weak coupling the gain in free energy from condensation grows faster with the magnetic field than 
this cost. As a consequence, at weak coupling a magnetic field can only break, never restore, chiral symmetry at a 
given chemical potential. This manifestation of MC is recovered in the Sakai-Sugimoto model for very large magnetic 
fields. For smaller magnetic fields, however, the opposite IMC is possible since the gain from condensation turns out 
to be less efficient in compensating the cost. This effect seems to be of general nature since it has also been observed 
for instance in NJL calculations (TJ [T3, HH, HI]. Nevertheless, our previous study was incomplete since we have for 
simplicity ignored baryonic matter. It did contain nonzero baryon number induced by a meson supercurrent via the 
axial anomaly (331436} but not "normal" baryons in the following sense. 

Baryons are introduced in the Sakai-Sugimoto model as D4-branes wrapped on the S 4 of the background geometry 
[l9l [33 , implementing the general idea of baryons in the AdS / CFT context [H, |39[ . Because of the flux carried by 
these D4-branes, each necessarily comes with N c string endpoints. The other end of each string is attached to the 
flavor D8-branes. In this way the AdS/CFT baryon realizes the usual picture of a baryon as a colorless object made of 
N c quarks which carry color and flavor. In the non-supersymmetric Sakai-Sugimoto model, the strings pull the baryon 
D4- branes towards and "into" the D8-branes [40] . Thus it turns out that a baryon can equivalently be described as a 
gauge field configuration with nontrivial topological charge in the world-volume gauge theory on the D8-branes. The 
charge of the wrapped D4-brane is related to the topological charge carried by the instanton [4l|. The latter, in turn, 
can be interpreted as baryon charge since it couples to the U(1)b part of the gauge group in the bulk, associated with 
the global group of baryon number conservation in the corresponding field theory at the boundary. The interpretation 
of baryons as instantons yields a natural connection to the Skyrme model where baryons appear as solitons |42j . (This 
connection is discussed in ref. [43| in the context of a five-dimensional gauge theory.) And indeed, it has been shown 
that the Sakai-Sugimoto model includes the Skyrme model in a certain limit [l9| . 

Spectra and properties of single baryons in the Sakai-Sugimoto model have been computed with a Yang-Mills action 
approximating the Dirac-Born-Infeld (DBI) action associated with the D8-branes 0, [37| (for analogous calculations 
in the non- antipodal version of the model see (44|). This approximation is not sufficient for our purpose since certain 
effects in a background magnetic field require the full DBI action, for instance the phase transition in chirally restored 
quark matter which is reminiscent of a transition into the lowest Landau level [27l . Eq . Using the DBI action, the 
description of baryons as instantons of a non-abelian gauge theory becomes extremely complicated. We thus restrict 
ourselves to single-flavor physics where the gauge theory in the bulk becomes abelian. [Strictly speaking, at Nf = 1 
there is no spontaneous chiral symmetry breaking because the U(1)a is always broken due to the axial anomaly; 
however, at large N c this effect of the anomaly is suppressed.] In this case, baryons have been approximated as 
pointlike instantons or, equivalently, pointlike D4-branes in the antipodal and non-antipodal versions of the Sakai- 
Sugimoto model, see refs. [46| and 1471 . respectively. In this sense one can consider the main part of the present study 
as an extension of the results of ref. [47| to nonzero magnetic fields. Baryons in a magnetic field have been considered 
previously j35j , however in the confined geometry with antipodal separation. Therefore, our study is also an extension 
of ref. [3a ] to the non-antipodal deconfined case, where the phase structure is much richer, as we shall see. Our study 
has some overlap with refs. (48l - [50| . While these references treat the baryonic phase purely numerically, we shall 
present analytical approximations which enable us to interpret the results physically. More importantly, however, we 
disagree with several results in these works, see footnote [5] in scc. HIIBl 
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C. Holographic vs. real-world baryonic matter 



Using the holographic setup just described to determine the effect of baryonic matter on the chiral phase transition 
and on IMC requires us to ask whether we expect to be able to draw any conclusions for QCD. It is known that 
baryonic matter in the Sakai-Sugimoto model is different from real-world baryonic matter in at least one important 
aspect. Namely, holographic baryonic matter is not self-bound, i.e., not stable at zero pressure. This unrealistic 
property originates from the baryon-baryon potential being repulsive for all distances [5l|. As a consequence, the 
zero-temperature onset of baryons is a second-order phase transition at the baryon mass hb = tub, and baryonic 
matter right after this onset has infinitesimally small density 0, In contrast, the onset of real- world baryonic 
matter is a first-order transition at fis = ms + Eq with the binding energy Eq ~ — 16MeV, and the density right 
after the onset is the nuclear ground state density no ~ 0.15 fm _ . 

In order to work out the differences between holographic vs. real- world baryonic matter in our context, we compute 
the baryon onset in the mean-field approximation of the Walccka model |52j in the presence of a magnetic field for 
electrically charged and neutral baryons and for various values of the anomalous magnetic moment. The Walecka 
model is a relativistic, field-theoretical model usually employed to describe dense nuclear matter in the context of 
neutron stars. It is constructed to reproduce properties of nuclear matter at the ground state density n . (To describe 
nuclear matter at densities beyond uq, a plethora of models exists, and it is an ongoing effort to exclude models for 
instance with the help of neutron star data.) In the simplest version of the Walecka model, baryon interactions are 
modelled by the exchange of the scalar sigma meson and the vector omega meson. Since the scalar meson exchange is 
responsible for the attractive interaction, the absence of a binding energy can be attributed to the absence (heaviness) 
of the sigma. And indeed, in the Sakai-Sugimoto model the lightest scalar meson is, in contrast to the real world, 
much heavier than the lightest vector meson. It is not completely understood whether this property is a large- N c 
effect or whether the strong coupling limit and/or a model artifact also play a role. It seems that the large N c limit is 
at least partially responsible since generalizations of the Walecka model to large values of N c have shown that nuclear 
matter tends to become unbound essentially as soon as N c is larger than three [lH, [HJ . 

Studies of nuclear matter in a magnetic field using Walccka-like models have been done previously in the context 
of neutron stars (55l - [57| . In these works, the conditions of beta-equilibrium and electric neutrality play a crucial role 
and thus the results are relatively complicated. Since we are interested in simpler, more fundamental questions we do 
not take into account these conditions. As a consequence, our results reflect the "pure" effect of a magnetic field on 
a single baryon species with a given electric charge and anomalous magnetic moment. They are thus, although less 
realistic, also interesting on their own beyond being a comparison to the holographic model. 

The remainder of the paper is organized as follows. In sec. |TT]we present the Sakai-Sugimoto framework and the 
solution of the equations of motion for the baryonic phase. These solutions are then used in sec. IHII to compute the 
transition between mesonic and baryonic phases in the plane of magnetic field and chemical potential. In particular, 
we shall discuss the role of the meson supercurrent for this transition. Sec. IIVI is devoted to the calculation of the 
onset of baryonic matter in a magnetic field within the Walccka model. In the final section of the main part, sec.|Vl we 
come back to our holographic approach and present the zero-temperature phase diagram, containing the chiral phase 
transition in the presence of baryonic matter. Readers who are only interested in the main results should consult fig. 
|5]and the brief summary of our observations in sec. IV Bl Sec. IVII summarizes all results and gives a brief outlook for 
future work. 



II. SETUP 



In this section we discuss one-flavor hadronic matter in a background magnetic field in the deconfined, chirally 
broken phase of the Sakai-Sugimoto model. The ingredients for our starting point can all be found in previous works 
wherefore we shall be very brief in the derivation of the action and the equations of motion. Details about the mesonic 
part in the same notation as used here can be found in ref. [13] while for the baryonic part we closely follow re f. [47| . 
For more general introductions into the Sakai-Sugimoto model see the original works 0, [2(| or the review [58l| or 
applications of the model for instance in refs. [59l-61|. 



A. Action 



The action we consider in the following is 

S = S'q + Sb ■ 



(1) 
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Here Sb is the contribution from the baryon D4-branes wrapped on the S 14 , and S' Q is the gauge field action on the 
D8-branes without the baryonic terms. Since we consider one flavor there is one D8-brane and one D8-brane. They 
correspond to left- and right-handed fcrmions, respectively. In the chirally broken phase considered in this section, the 
two branes are connected and have a nontrivial embedding in the background geometry, which has to be determined 
dynamically. The prime denotes a modification to the Chern-Simons (CS) part of the action by certain boundary 
terms. These terms have been argued to be necessary in the presence of a magnetic field [35| and have been discussed 
further in ref. (62|, where in particular their problems in the context of the chiral magnetic effect have been pointed 
out. Starting without this modification we have 



So — Sdbi + Sq S . 



(2) 



The notation S^g indicates that there is a second CS contribution from the baryons which we shall call Sq§ , see eq. 
(fTTjl . The DBI part is 



drd 3 x I dC/e-* v /det(,9 + 27ra'F). 



(3) 



9 is the D8-brane tension, V± = 8?r 2 /3 the S 4 volume, and e* = g s {U/Rf/ 4 
and the curvature radius R of the background geometry. The background 



Here g is the induced metric on the D8-branes in the deconfincd geometry (see eq. ©), F the field strength tensor of 
the U(Nf) gauge theory on the branes, which in our case is abelian since we work at Nf = 1. Furthermore, a' = £ 2 
with the string length £ s , T s = (27r)~ 8 < 
the dilaton with the string coupling g 

geometry is determined by N c D4-branes (not to be confused with the baryon D4-branes wrapping the S 4 ) and 
fixed throughout the paper since we work in the probe brane approximation. The integration is over Euclidean 
four-dimensional space-time with imaginary time r and over one half of the connected flavor branes parametrized by 
the holographic coordinate U <E [U c ,oo], where U c is the tip of the U-shaped, connected branes and U = oo is the 
holographic boundary where the gauge symmetry becomes global and corresponds to the chiral symmetry group. In 
the presence of pointlike baryons, as discussed below, it turns out that there is a cusp in the profile of the branes at 
the tip U = U c . The CS part in the gauge A\j = is 



c(0) 



24tt 2 



dr d A 



dU A^FuvFpvt 



(4) 



where jx, v, . . . = 0, 1, 2, 3 and 60123 = +1. 

Our ansatz for the gauge fields takes into account the nonzero boundary values Aq(U = 00) = fi q and A±(U — 
00) = —xiB where fi q is the quark chemical potential and B the modulus of the non-dynamical background magnetic 
field, pointing in the 3-direction. This ansatz requires a nonzero A 3 (U) to fulfil the equations of motion, and even a 
nonzero boundary value A^(U = 00) to minimize the free ener gy. This boundary value is interpreted as the gradient of 
a meson field which corresponds to a meson supercurrent [3J-|36[ . Within this ansatz and including the modification 
of the CS action mentioned above, we exactly follow the steps in ref. [13] to obtain 



2 T 



du Cq , 



(5) 



where V is the three- volume, T the temperature, and we have abbreviated 

N c R 2 



Af 



6tt 2 (27ra') 



A3 



and 



(6) 



Co 



y/(u* + & 2 « 2 )(1 + fa' 2 - a' 2 + vFfxf) + - b(a 3 a' - aoa' 3 ) . 



(7) 



with the first and second term coming from the DBI and CS contribution, respectively. We have introduced the 



following dimcnsionless quantities: the gauge fields a M = ^j^-A^, the magnetic field b = 2-na'B, the holographic 
coordinate u = ^, and the coordinate of the compactified extra dimension X4, = X± € [0, ^ ] with the Kaluza- 

Klein mass Mkk- The D8- and D8-branes are asymptotically, i.e., for u — > 00, separated in the X4 direction by a 
distance L which we shall take to be small with respect to the maximal separation on antipodal points of the circle, 



L <C m ■ Moreover, the prime denotes derivative with respect to u, and 

47T. 



/(«) 



f 



If 



a 



1/2 



-TR. 



(8) 
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To obtain eq. ([5]) we have used the relation R 3 = irg s N c £ 3 and the induced metric on the D8-brancs 



ds 2 = u 3 / 2 [f(u)dT 2 +S tj dx l dx3}+ R 



u 3 / 2 f(u) 



u 3 ' 2 x' 2 



du 2 + u l l 2 dVL\ \ , (9) 



S D4 + 4s ■ ( 10 ) 



where dfl 2 is the metric of the S" 4 . 
Now we turn to the baryonic part, 



To understand the origin of S^s we mus t take a detour via the nonabelian case Nf > 1. In that case, eq. ([!]) is 
not the only CS contribution. The integrand of eq. ^ is of the form A F 2 with a gauge field A and field strength 
F from the U(l) part of the bulk gauge group U(Nf) = SU(Nf) x U(l). In general, there is also a contribution 
where the integrand has the form ATr F 2 where A is from the abelian part, but F from the nonabelian part with the 
trace taken over flavor space (see for instance Refs. [36|, H3] for an explicit form of the CS action in terms of abelian 
and nonabelian parts). Baryon number is generated from the term which couples to A$, i.e., we are interested in a 
contribution of the form 

8& = £[ Ao Tr F 2 . (11) 



xU 



Here the integration goes over Euclidean space-time R and the holographic direction U, parametrized by 
(t, x , x 2 , x 3 , u). From this expression we see that the instanton charge, 



Na = ~^\ TrF 2 , (12) 



8- 



R 3 xW 



with integration over three-dimensional space and holographic direction, parametrized by (x 1 , x 2 , x 3 , u), is identical 
to the baryon charge. For Nf = 1, there is no non-singular instanton solution, and we use a pointlikc approximation 
localized at u = and homogeneous in three-space p7| (for attempts towards inhomogeneous solutions see for 
instance refs. [63l ItMj). 



Thus the CS contribution becomes 



-LtyF 2 = -^5{u-u c )d 3 xdu. (13) 
on V 



S { cl = ~Mn iaQ {u c ), (14) 



where we have defined the dimcnsionless baryon charge density 

N C N 4 R 

n4= AT^w- (15) 

The first term on the right-hand side of eq. (fTU|) is the action of \N&\ D4-branes wrapped on the S 4 , corresponding 
to |i\?4 1 (anti-)baryons. In the case of anti-baryons, these are D4-branes with negative charge N4 < 0. In our 
approximation of pointlike D4-branes, 

S D 4 = \N 4 \T 4 J dn 4 dTe-*^U/ = ^AA^V/K) l"4| , (16) 

where T4 = (27r) _4 £j 5 is the D4-brane tension and where, in the second step, we have used the explicit forms of the 
metric and the dilaton. For later use it is convenient to write Sd4 in an alternative way, introducing the (dimcnsionful) 
mass M q of a constituent quark within the baryon, 

ci_jii M M ^ c 

TV q q 2na' 3 

such that N c M q is the baryon mass. 

We can now put together the contributions ([5]), (fl4|). and (p~6]> to obtain the total action 

V f°° V 



Sd4 = ^^^M 9 , M^^^^TK), (17) 



T J„_ u T 

with Cq given in eq. ([7]). 



S=Af— I du£ Q +Af—?i 4 sgn(n 4 )y y/ f{u c ) - a„(u c ) , (18) 
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B. Equations of motion 



From now on we shall set f(u) ~ 1. This simplifies the analysis of the chirally broken phase tremendously since only 
within this approximation we can find semi-analytical solutions. As a consequence, our results are more transparent 
compared to a purely numerical analysis and the physical interpretation becomes simpler. This approximation can 
be read in two ways because f(u) ~ 1 is equivalent to ut <C u c (since u > u c ), and ut and u c are determined by 
independent scales of the model. From eq. © we see that 117- is proportional to the temperature. Hence f(u) ~ 1 can 
be interpreted as the zero-temperature limit. We need to remember, however, that in the Sakai-Sugimoto model the 
confined phase is the ground state for T <T C = . Therefore, we can only take the T — > limit of the deconfincd 
phase if simultaneously we also let Mkk — > 0, i.e., let the radius of the compactified extra dimension go to infinity. 
The second option to fulfil ut <C u c is to change u c at fixed Ut- This can be achieved by decreasing the asymptotic 
separation of the D8- and D8-branes, i.e., by letting L — > 0. Geometrically, however, this is not so different from the 
first limit: increasing the radius of the circle along which the branes are separated while keeping their distance fixed 
corresponds to decreasing the separation while keeping the radius fixed. In either case we arrive at a "decompactificd" 
limit where the gluon dynamics is decoupled from the system and which is not unlike a nonlocal NJL model on the 
field theory side [H, [H] . 

The following parts, including the calculation in appendix [3 and the next subsection III C[ are the most technical 
ones of this paper. They serve to collect all necessary results which shall be evaluated and discussed in sees. IHII and 
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With f(u) ~ 1 the equations of motion become in the integrated form 



y/\ + a'i - a' 2 + u 3 xf 



y/1 + - a' 2 + u 3 xf 



,a/u 5 + i 



y/1 + a 12 - a'J + uH'i 



= 3603 + ri4 , (19a) 
= 3ba + d, (19b) 
= k. (19c) 



Here, d and k are integration constants, to be determined. The integration constant of the first equation is fixed to 
be ri4 because of the pointlikc charge sitting at u = u c . As written, these equations are identical to the ones without 
baryonic matter. The singularity at the tip of the brane, however, gives rise to different boundary conditions. 

Without explicitly solving the equations of motion we may compute the nonzero components of the four-current 
jv- . The current is given by 



J 1 



dC 



(20) 



since the baryonic part Sb does not contribute. (Here we have used the dimensionful gauge field to obtain a 
dimensionful current.) This yields 

J° = (ffc + n 4 ) , J" = - 2 -^N fa + d) , (21) 

where we have used the equations of motion and where fi = do (00) is the dimcnsionless quark chemical potential, 
ji = Mq; an d 3 = 03(00) is the dimensionless meson supercurrent which has to be determined dynamically from 
minimization of the free energy. The dimensionless baryon density, 

3 

n=-b] + n Al (22) 

thus receives contributions from the meson supercurrent and from "normal" baryons . In our model, both contribu- 
tions are of topological nature since they both originate from the CS term. We shall see below that n is indeed the 



1 Note that the factor 27 ^* Af, which turns n into a dimensionful density, scales with N c (while n does not). Therefore, although it may 
seem a bit confusing, n is a (dimensionless) baryon density, while S° is the (dimensionful) quark density. (In our previous work [27| the 
terminology was more sloppy since we referred to n as a dimensionless quark density.) 
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baryon number density also in the thcrmodynamical sense, obtained from taking the derivative of the free energy with 
respect to /x, see eq. f|43|) . This consistency is, in the presence of a magnetic field, only achieved by the modification 
So —> S' introduced in ref. [35|. Analogously, J 3 can also be obtained by taking the derivative of the free energy 
with respect to the corresponding source, which in this case is the supercurrent j. Consequently, the condition that 
j minimize the free energy is equivalent to J 3 = 0. Nevertheless, one should not think of J 3 = as an externally 
imposed constraint; the system rather finds its minimal free energy for any given /i and b, and J 3 — is a mere 
consequence. Using J 3 = in Eq. (|21[) , we can immediately read off 

d = -\bfi. (23) 

We also need to minimize the on-shell action with respect to the parameters 714 and u c . By using partial integration 
and the equations of motion, we have 

dC ddj | dC dx 4 j ± u c ^ ^ da (u c ) 

^ da'; dn4 Dx'a dn 4 I 3 c dn 4 

i=0,3 1 * / 

/ u—u c 

where the upper (lower) sign holds for 7i 4 > (n 4 < 0). The derivatives involving X4 vanish since we keep the 
asymptotic separation of the flavor branes fixed, 2:4(00) — X4(u c ) = | with the dimensionless separation I = The 
derivatives involving the gauge fields vanish at u = 00 since ao(oo) = /i, 03(00) = j are held fixed. We also impose 
&3(uc) = such that the only nonvanishing term from the parentheses comes from the derivative involving ao, at 
u = u c . The equation then reduces to the simple condition 

ao(ue) = ±y. (25) 

The minimization with respect to u c becomes 

where we have used eq. (|25j). Here we have to be more careful at the lower boundary because at this boundary ao, 03 
and 24 depend on u c explicitly as well as through the variable u. Therefore, 



da Q (u c ) , da (u) 
- a (u c ) + 



8u r 3u r 



(27) 



and the same for 0,3 and X4. Here, the left-hand side denotes the total derivative with respect to u c . In the first 
term on the right-hand side the derivative acts on the dependence through u, while in the second term it acts on the 
explicit dependence. Making use of this relation and of the equations of motion, eq. (ffij)) can after some algebra be 
written as 

Kl VffK) , . 

-ir = ^' (28) 

where we have abbreviated 

g{u) = u s + b 2 u 5 - k 2 - u 3 [(bu c ± d) 2 - nl] . (29) 



The condition (|28j) can be interpreted as a force balance equation which says that the force from the D8-brane tension 
cancels the one from the weight of the baryon D4-branes [47| . This equation is now used to determine k, 

^ = „s + ^_„. [(J „ c±(i ) 2 - n a-„3(^. (30) 

In appendix [S] we discuss the explicit solution of the equations of motion in detail. Here we only summarize the 
results. The gauge fields become 

. fi / ^ sinhy \ ^ //i u c \ ^ / sinhj/ coshy \ . 
2 \ sinhi/oo/ \2 3/ 00 ysinhyoo coshj/oo/ 

^coshy-1 fl-iu c \ /coshy-1 sinhy \ 

a%{y) = o^-r + ( o =F -5-] coshj/oo — — , (31b) 

2 sinn yoo \ 2 3 / \ Sinn j/oo cosn yao J 
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with the new variable 

y(u) = 3b 

and i/oo = y{u = oo). The embedding of the flavor branes is 

xa{u) = k 

With d from eq. (j2"3")l and 



V dv 



dv 



v 3/2 ^/g{v) 



3b 



714 



2 sinh y a 



/' 



2u r 



cosh y 



(32) 



(33) 



(34) 



k and thus g{u) can be expressed in terms of u c and yoo which have to be determined numerically from the two 
coupled equations 



ih 



3b 



00 u^ 2 du 



du 



(35) 



Since yoo is manifestly positive, eq. (f3"4l shows that /i and 714 must have the same sign, as it should be. Note also that 
in k and g(u) the sign of 714 appears only in a product with fi. Therefore, the equations for y^ and u c only depend 
on We may thus restrict ourselves in the following to positive n 4 without loss of generality. 
Finally, the supercurrent becomes 



u c cosh y a 



while the total density is 



n = 3b I u — 



3 sinh yoo 

u c \ 1 + cosh y c 



2 sinhy 



(36) 



(37) 



C. Free energy 

The free energy of the baryonic phase is 

rp POO 

fiv = -^Son-shcii = W / du£ , (38) 

since the baryonic part Sb does not contribute, due to eq. (|25|) . (Of course, the effect of the baryonic part is implicitly 
present through the equations of motion.) By inserting the explicit solution - most conveniently in the form of eqs. 
HI) - we find 

n v = nT duu^ uh + bV ~ ^ + d? ~ ni] +N (™ - ^) . (39) 

This form of the free energy is useful for a comparison with the purely mesonic phase, see eq. (3.14) in ref. [27j |. For 
a proof that the derivative of f2y with respect to /1 is indeed the (negative of the) total density we rewrite 



f°° , \J~gJu) . k£ yoo 2 2 fin n 4 u c 

jr = J Uc du ^ + T + ^ [(5uc+d) -" 4] "T + ^- (40) 

Now we note that f2 v = Qv[fJ-, u c (fi, 6), y o{fJ', b)] and take the derivative with respect to fi at fixed 6, including the 
implicit dependence through u c and yoo. We obtain 

1 dfl\/ ri4 du c (6it c + d) 2 — n\ dyoo , d ( fin | n^u c \ 



J\f dfi 3 dfi 66 dfi dfi 
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FIG. 1: Solid lines: baryon density n.4 (left) and location of the tip of the connected flavor branes u c (right) at vanishing 
magnetic field. The baryonic phase becomes favored for /i > u c /3, indicated by the dotted line in the right panel. For chemical 
potentials smaller than that value, the mesonic phase is the ground state, in which - at vanishing magnetic field - the baryon 
density vanishes and u c is constant. The dotted line in the left panel is the linear approximation (|47(l . Here and in all other 
plots we have set I = 1 for notational convenience. The dependence on £ is always simple and corresponds to a rescaling of the 
axes, in this plot — > l 2 ^j,, u c —t l 2 u c , m — > £ 5 n4- 



With the explicit expressions for Ua and n in eqs. (|34|) and (|37l) one computes 



d f fin n i u c \ riiduc (bu c + d) 2 - n\ dy^ 

n +^TlTT SI oTT ' ( 42 ) 



dfj, V 2 6 / 3 dfi 66 dfi 

such that indeed 

dn v 2na' 



dfXq R 



Nn= J°. (43) 



III. HOLOGRAPHIC BARYON ONSET 



A. Zero magnetic field 

Let us start with the simplest case of a vanishing magnetic field. The numerical solution of eqs. ([33]) for this case 
yields the results shown in fig. [1] For chemical potentials smaller than ^ no solution with n 4 ^ exists, and the 
mesonic phase is the ground state. This phase has already been discussed in ref. (27j and is obtained by omitting 
the baryonic contribution Sb in the action ([1]). To find the transition between the mesonic and baryonic phases, the 
"baryon onset" , we consider the respective solutions wherever they exist and compare the respective free energies. 
(The mesonic and baryonic phases are not distinguished by any symmetry, the baryon onset is only a phase transition 
in the sense that the second derivative of the free energy is discontinuous.) The numerical calculation shows that as 
soon as a baryonic phase with 72,4 > exists it has lower free energy than the mesonic phase. This statement holds 
for all magnetic fields b. 

Before we come to the onset at finite b in detail, we discuss the results for the 6 = case analytically. From the 
numerical solution we read off that yoo is linear in 6 for small 6. Thus we insert the ansatz = 36?^, with to 
be determined, into eqs. (|3~i)) and (f3"T)l to find 

n 4 (6 = 0) = n(b = 0) = ^— =— ^ , (44) 

where = 27 ^* M® is the dimensionless constituent quark mass at vanishing magnetic field, mP q — ^ with u a c = 
u c (b = 0). Note that yj^, u° c and thus also m° are still complicated functions of fx, to be computed numerically in 
general. From eq. (HU) we conclude that a baryon chemical potential at least as large as the baryon mass is needed 
to obtain a nonzero baryon density. This shows that there is no binding energy for our holographic baryonic matter, 
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quark density N c N/t/V 


B = 


2 A £ 


B — > oo 





TABLE I: Quark density N c Ni/V induced by wrapped D4-branes close to the baryon onset at vanishing and at asymptotically 
large magnetic field B. Both expressions depend linearly on the difference between the quark chemical potential and the 
(effective) constituent quark mass, indicating a bosonic nature of the holographic large-iV c baryons discussed here. Since, at 
nonzero magnetic field, the effective constituent quark mass in a baryon is modified by the meson supercurrent, the baryon 
onset for B — ¥ oo is at 2M^° , not at . For the values of the numerical constants £, £ see eqs. (|48j) and (1531) . 



since then the onset chemical potential would be lower than the mass, as it is the case for real- world nuclear matter. 
As a consequence, our baryon onset is a second-order, not first-order, phase transition. We shall see that this is also 
true for nonvanishing magnetic fields. 

In the vicinity of the onset we may compute the density analytically. First, by inserting n 4 — into eqs. (|35p . we 
obtain at the onset, 



c, onset 



16tt 



. 16/ 



. 16/ 



e r 



5127rr 



. 16/ 



'IT 

. 16/ 



r(&) 
r(&) 



(45) 



and thus 



r-onset " l q 



0.17 



(46) 



This result is in agreement with rcf. [42j, see for instance fig. 10 in this reference. For chemical potentials larger than 
but close to fi® nset , the dimensionless density behaves as (see appendix IB II for the derivation) 



n 4 {b = 0) 



(/i-m°) + 0[(M-m°) 2 ], 



(47) 



where 



r(A)r' 



16/ 



r(£)r(g) 



o.n 



(48) 



We plot the linear approximation (|4T|) in comparison to the full result in the left panel of fig. [TJ 

This linear behavior might seem unexpected because for fermions in 3+1 dimensions, at least in the non-interacting 
case, one would have rz.4 ~ (fx — to) 3 / 2 just above the onset. It is therefore interesting to compare 714 from eq. (|4"7| 
with the density of a Bose condensate instead. To this end, consider a 4 model with chemical potential at zero 
temperature, f2 = — ~ M <f) 2 + j<ft 4 - Minimizing f2 with respect to the condensate <fi and inserting the result into the 
density n = yields n = /i M = — to) + — to) 2 ]. The linear term looks exactly like the one in eq. 

(|4T|) (in the <p A model, the quadratic term is positive; this is in contrast to our holographic result, as fig. Q] shows). 
The similarity of the baryon density with the density of a Bose condensate is not too surprising because we work in 
the large- A^ c limit. In this limit, changing the number of constituent quarks from even to odd, i.e., from A^ c to N c + 1 
does not make a difference. Therefore, in our context one may very well talk about a condensate of baryons. This 
terminology is for instance used in ref. [46j. It is instructive to express eq. ()47|) in terms of dimensionful quantities, see 
first row of table HI where we show the dimensionful quark density J® . In the case B = this density only receives a 
contribution from "normal" baryons, J°{B = 0) = N c Ni/V, see eqs. (|15p and (p?Tj) . Comparing this expression with 
the </> 4 model, we may identify M^ ZeZMss with a coupling constant for an effective repulsive interaction between the 
baryons near the onset. We see that this coupling is proportional to the 't Hooft coupling A. Interestingly, it gets 
strong for small asymptotic separations L of the flavor branes, measured relative to the maximal separation tt/Mkk- 
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B. Nonzero magnetic field 

After having checked numerically that baryonic matter, where it exists, has always lower free energy than the 
mesonic phase, we know that the baryon onset line in the diagram is defined by n 4 = 0. With the help of eq. 
()34p this gives a condition for coshy,^, into which we insert y^ from eq. (|35[) . Together with the condition for the 
asymptotic separation from eq. (|35|) this yields two coupled equations for u and u c , 

3/i , [°° Qbu c u 3 ' 2 du 
= cosh / — ^^^^^^^^^^^^^=^^^^^^^^^^^^^= , (49a) 

2Mc ~ 3i " Jl yj 4u 2 [(u s - l)u 3 c + (u 5 - 1)6 2 ] + {u 3 - l)b 2 (2u c - 3uf 



u 



1/2 1 f °° V 4w c(«? + b 2 ) - b 2 (2u c ~ 3/i) 2 du 



(49b) 



u 3 / 2 y / 4u2[( w 8 - i) u 3 + (u 5 - 1)6 2 ] + (u 3 - l)b 2 (2u c - 3/i) 2 ' 
where we have changed the integration variable u — > -7- . Since the right-hand side of eq. (|49aj) is > 1 we conclude 

y < /Wet < — ■ (50) 

(Remember that u c itself depends on b and /i.) We have seen above that the lower boundary is saturated for 6 = 0. 
The upper boundary is saturated for b — > oo, as can be seen from setting 77,4 = in eq. (|34p and using — > 00 at 
asymptotically large b. Hence, by inserting /x = -j^ into eq. (|49b[) we find the asymptotic u C! which in turn gives the 
asymptotic onset 



32tt 
3^~ 



Closely above that value the behavior of the asymptotic density of "normal" baryons 114 is given by (for a derivation 
see appendix IB 2[) 

n 4 (b = co) = y / , g +0[( M -2m g °°) 2 ], (52) 

with 

2 r(J-) 

f = l ^—^ff~ 0.5194. (53) 



150F r 



The corresponding dimensionful quark density is shown in the second row of table [U (Remember that the total baryon 
density n also receives a contribution from the supercurrent, n = |6j + 77,4 with j(b = 00) = u c /3.) Interestingly, at 
asymptotically large magnetic fields the explicit dependence on the model specific constants L, A and A/kk drops out. 



This kind of "universality" has already been observed in the phases without baryons, see for instance eq. (4.9) in [27 1. 

One might have thought that the baryon onset always happens at fj, = ^ since N c ^j- is the baryon mass given 
by the action of the D4-branes. To understand why the result deviates from this expectation, it is instructive to 
consider the (unphysical) case of an isotropic meson condensate where the meson supercurrent vanishes, j = 0. In 
this "cleaner" case, the true vacuum with zero pressure P = and zero baryon density n = is the ground state 
below the baryon onset. Now, cq. (|A9b[) shows that the onset indeed occurs at n = ^ for arbitrary values of the 
magnetic field. As a consequence, the chemical potential equals the energy per baryon e/n at the onset (this follows 
from the thermodynamic relation P = — e + (in and P = at the onset). The situation with a nonzero supercurrent 
j is different. Here the pressure and the baryon density become nonzero as soon as we switch on /i, and n is always 
larger than e/n. Then, at some point, \i is large enough to support the baryon number induced by j and by "normal" 
baryons. This costs more energy than having only "normal" baryons and thus the onset happens later than for j — 0. 
We compare both onset lines, the unphysical one from j = (dashed) and the physical one with j ^ minimizing the 
free energy (solid), in fig. [5] 

The difference between the two lines can also be thought of as follows. Between the dashed and solid lines the system 
might "think about" adding baryons because there is enough energy available from the chemical potential. If the 
system decided to do so, it would at the same time have to decrease the supercurrent discontinuously to have enough 
energy available for the baryons. (The dashed line is the extreme case where it would have to force the supercurrent 
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2.0 
1.5 
£> 1.0 
0.5 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 2: Onset of baryonic matter (solid line) in the plane of quark chemical potential ^ and magnetic field 6. For comparison, 
the dashed line indicates the baryon onset in the unphysical case of a vanishing meson supercurrent. In this case, the onset 
occurs at ^ for all 6. Due to the meson supercurrent, the effective baryon mass is more complicated and the onset is somewhere 
between ^ (b = 0) and 2|£ (b = oo). The arrow indicates the asymptotic value of the onset line given in eq. (|51|l . In the 
present model, the baryon onset is always a second-order phase transition, in contrast to real nuclear matter. For now, we have 
ignored the chirally restored phase. We shall see in sec. [V] that for sufficiently large magnetic field, b > 0.25, the transition to 
baryonic matter is replaced by a transition to chirally symmetric quark matter. 



to zero.) This, however, would lead to a decrease in the total baryon density upon increasing the chemical potential, 
which is a thermodynamically unstable situation. (At the dashed line, the baryon density would jump to zero because 
the energy is just enough to start adding baryons 2 .) For any given b the solid line marks the "earliest" possible point 
where baryons can be put into the system without such a thermodynamic instability. 

What can we learn about real- world baryonic matter from these results? We already know that the nature of the 
onset is different due to the lack of binding energy. But why do our holographic baryons effectively become heavier 
in a magnetic field? Can we draw any conclusions from this observation? Within our model, the increasing critical 
chemical potential has two reasons. First, as just discussed, it is the meson gradient which is responsible for increasing 
the baryon onset from ^ to Second , Ur itself increases with the magnetic field. In the mesonic phase u c can be 
interpreted as a constituent quark mass [60, 65]. In the chiral limit of vanishing bare quark masses, it is the chiral 
condensate which induces such a mass. Therefore, the increase of u c in a magnetic field suggests an increase of the 
chiral condensate and thus is a manifestation of magnetic catalysis (MC). The constituent quark mass in a baryon is 
different from the constituent mass in the mesonic phase, ^ vs. u c , but both are proportional to u c . Consequently, 
MC seems to be responsible for the heaviness of magnetized baryons. 

IV. BARYON ONSET IN A RELATIVISTIC MEAN-FIELD MODEL 

In this section we employ the Walecka model ,.52j at zero temperature in a background magnetic field. The Walecka 
model is a relativistic model for dense nuclear matter, where nucleons (or, in extensions of the model, hyperons) 
interact via Yukawa exchange of mesons. In the simplest, isospin-symmetric, version considered here, nucleons interact 
through the scalar sigma meson and the vector omega meson. This is sufficient to model the realistic nucleon-nucleon 
interaction which is known to have a repulsive short-range (omega) and an attractive intermediate and long-range 
(sigma) part. Nuclear matter being stable at zero pressure, having a finite binding energy Eq ~ — 16.3MeV, a 
saturation density no — 0.153 fm -3 , and showing a first-order liquid-gas phase transition are all manifestations of 
this simple but crucial property of the interaction. Beyond the saturation density, the properties of nuclear matter 
are poorly known, and the Walecka model is only one of many competing models. However, for our comparison to 
holographic baryonic matter, we are primarily interested in the baryon onset, and we know that for this purpose - at 
least for vanishing magnetic field the Walecka model describes, by construction, real-world nuclear matter. 



2 This unphysical onset is considered in ref. l50j : the T = onset in fig. 9 of this reference corresponds to the dashed line in our fig. [2] 
It seems that this discrepancy originates from incorrectly scaling 5d4 with the total density n, not 714, see eq. (2.4) in ref. [4H which 
apparently is used in ref. |5C| (and in [4H I. This difference is very important for the topology of the resulting phase diagram: as we 
shall see in sec. IV Al the physical onset line intersects the chiral phase transition line, see fig. [6] this is not the case for the dashed line. 
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A. Lagrangian 



The following setup is essentially taken from works where dense nuclear matter in a magnetic field has been 
considered in the astrophysical context, see refs. [55l - l57| and references therein. In neutron stars, the simplest version 
of baryonic matter consists of neutrons, protons, and electrons. Their various chemical potentials are related through 
the conditions of beta equilibrium and (global) electric charge neutrality. For our purpose, these complications are 
irrelevant since for a comparison to our holographic results we are interested in the behavior of a single baryon species 
with a given electric charge and a single baryon chemical potential in a background magnetic field. Therefore, our 
results will be simpler and more transparent, albeit less realistic, compared to the results for astrophysical nuclear 
matter. 

We start from the Lagrangian 

C = C B + Ci + C M , (54) 
containing a baryonic part Cb, an interaction part £/, and a mesonic part Cm- The baryonic part is 

C B = 4> \vfDp - m B + MB7° - ^Kd^F^j V , (55) 

with the baryon mass raj and the baryon chemical potential /jb- The baryons feel the magnetic field through the 
covariant derivative = <9 M + iqA^ with the baryon electric charge q and the electromagnetic gauge field A^, 
which encodes the background magnetic field in the ^-direction, A^ = (0, X2B, 0, 0). In addition, baryons have an 
anomalous magnetic moment n whose effect is included by a magnetic dipole term with cr M „ = i [7^ , 7„] and the field 
strength tensor F^ v = <9 M A v — d v A^ L . It is important to keep in mind that this term is obviously an effective, not a 
fundamental, way to take into account the anomalous magnetic moment. In particular, the present approach can not 
be trusted for arbitrarily large magnetic fields, as we shall see more explicitly below. 
The interaction term consists of two Yukawa contributions for the a and the to, 

Ci = gai>o"ip - g^-i^u^ip , (56) 
with coupling constants g ai guj > 0. The mesonic part includes cubic and quartic scalar self-interactions, 

Cm = \(d^a - mla 2 ) - ifi"^ + \ml^ - h -m B {g a <jf - C -{g a af , (57) 

with f2 Ml/ = d^<jj v — d v Lo^ and the sigma and omega masses m a and . With the given values tub — 939 McV, 
m u = 783 MeV, m a ~ 550 McV, the model has four free parameters which are fitted to reproduce the saturation 

density, the binding energy, the compressibility, and the Landau mass at saturation (all for B = 0), which gives 

22 

I 2 - = 6.003, I s = 5.948, b = 7.950 x 10~ 3 , and c = 6.952 x 10~ 4 . We shall employ the mean-field approximation where 
the meson fluctuations are neglected, and the meson condensates a and uiq have to be determined from minimization 
of the free energy. The basic equations of the model in this approximation are as follows (see for instance ref. [66j for 
a pedagogical derivation). The pressure is 

P = -\™l° 2 - ^rn B {g^) 3 - \{g.of + + P B , (58) 

where, at zero magnetic field, the renormalized baryonic pressure is given by 



/ In 



l + e -(**-M.)/T 



(59) 



(2tt) 3 

Here, /i* = [Ib — .9^0 plays the role of the Fermi energy at zero temperature (the chemical potential in the ther- 
modynamic sense it still /is, not /x*). The baryon dispersion is e& = \/fc 2 + ml with an effective baryon mass 
to* = ttib — go®- The stationarity equations for the meson condensates are 

g 2 g 2 

m* = ni B — 4ri 8 + — v \bmB(mB — to*) 2 + c(tob — to*) 3 ] , (60a) 

TO J. TOJ. L J 

wo = -^n B , (60b) 

TO* 

where ns = and n s = — |^ are the baryon and scalar densities, respectively. 

Next we include the magnetic field, wherefore we need to distinguish between charged and neutral baryons. For 
our purpose the only relevant effect of the magnetic field concerns the single-baryon dispersion relations, which can 
be found in ref. [53 ■ 
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B. Charged baryons 

Our reference example of charged baryonic matter is pure proton matter with charge q = +e and an anomalous 
magnetic moment k = 1.79 /ijy where /U;v = 2^7- — 3.15 x 10~ 18 McV G _1 is the nuclear magneton. We shall, however, 
vary the anomalous magnetic moment below to study its effect on the baryon onset. Remember that we ignore any 
neutrality constraint and Coulomb effects. For charged baryons the dispersion has to be replaced by 



kf, + Ml s , M V)S = ^ml + 2v\q\B - snB , (61) 



with fc|| being the momentum in the direction of the magnetic field, s = ±1, and v = n + (1 — ssgnq)/2 with 
n = 0, 1, 2, . . ., such that for positive charge we have v = 0, 1, 2, . . . for s = +1 and v = 1, 2, 3, . . . for s = —1, and 
vice versa for negative charge. This means that the lowest Landau level (LLL) v = is only populated by s = +1 
fermions for q > and s = — 1 fermions for q < 0, while both s = +1 and s = — 1 contribute to all higher Landau 
levels. The three-dimensional momentum integral has to be replaced by a one-dimensional integral over fc|| and a sum 
over Landau levels and spin degrees of freedom, 



(2tt) 3 2tt 2 



s=± 







With these replacements in the baryonic pressure (|59[) we obtain at zero temperature 

\q\B 

S — ± V 



?b = PE - < s m ) , (63) 



and 



II B 



^ yy - ^ ■■■ / ' r - ' "• (64b) 



where kp tVtS = \/ fJ-'i — M% s is the longitudinal Fermi momentum, and where the Landau levels are occupied up to 



(11,+shB)^ — m, ^ inserting eqs. (|64[) into eqs. (|60p and solving the resulting equations for to*, cjo ; one can 

compute the thermodynamic properties for arbitrary (Ib and B. Since we are only interested in the baryon onset, we 
have the additional condition P = 0, where P is obtained by inserting eq. (|63| into eq. (|58|) . Hence we have a system 
of three equations to be solved for to*, wq, and at any given £>. 

The result for pure proton matter and three more (unphysical) types of charged baryonic matter with q = +e, 
distinguished by different values for k, is shown in fig. El where we plot the onset lines and, for proton matter, the 
baryon density along the onset. (In order to reproduce the nuclear ground state density n at B = we have multiplied 
the pressure by a factor of 2. In other words, we have started from the isospin-symmctric, B = Walccka model for 
protons and neutrons and have added the effect of the magnetic field as if both nuclear species had the same charge 
and anomalous magnetic moment.) 

The B = onset occurs at /is = ms + Eq ~ 922.7MeV where Eq ~ — 16.3MeV is the nuclear binding energy. 
Magnetic fields of the order of 10 18 G and larger change this significantly, in accordance with the simple estimate in 
sec. II At We observe oscillations due to the Landau level structure in the onset line as well as in the onset density. 
At sufficiently large magnetic fields only the LLL is occupied. (If there was no binding energy, the density at the 
onset would be infinitesimally small, and all along the onset line the LLL would be the only relevant state. In other 
words, only due to the presence of a finite binding energy the onset line shows a behavior reminiscent of de Haas- van 
Alphen oscillations at small magnetic fields. 3 ) Since we have, without loss of generality, fixed the electric charge to 



3 The higher Landau levels in fact induce cusps in some of the onset lines, see for instance the case n = where fJ. C g Bp — 920 MeV. In 
a small vicinity around this cusp there are, for a given B, two solutions for the onset: one where only the LLL becomes occupied, and 
one where the lowest and the first Landau level become occupied simultaneously. The resulting two onset lines intersect, and for any 
given B the line where the onset occurs at a lower fis is the physical one. The reason is that the "later" onset would keep the system 
in the vacuum state P = 0, in a region where the "earlier" onset already leads to P > 0. (For (j,b < f^ C g Sp , the LLL onset is "earlier", 
for fig > ffg^ it is "later", hence the cusp.) 
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Hb [MeV] B [10 iy G] 

FIG. 3: Left panel: zero-temperature transition in the plane of baryon chemical potential /is and magnetic field B from the 
vacuum (to the left of each line) to baryonic matter (to the right of each line) for charged baryons with q = +e and various 
values of the anomalous magnetic moment k (in units of the nuclear magneton pi at). Each line represents a first-order phase 
transition since the baryon density jumps from zero to a finite value no- Right panel: corresponding baryon density no along 
the onset line on a doubly logarithmic plot for the case of pure proton matter, k — +1.79 (the other three cases from the left 
plot would lead to similar curves). The oscillations are due to successive occupation of Landau levels. For B > 0.5 x 10 19 G 
only the LLL is occupied. 



be positive, the LLL is occupied by ,s = +1 baryons in all four cases considered here. Note that in one of the four 
cases shown in the figure, the anomalous magnetic moment n is negative. Since k < favors s = — 1 baryons, it 
can in principle be less costly to put baryons in the v = 1, s = — 1 state than in the v = 0, s = +1 state such that 
for sufficiently large magnetic fields all baryons sit in the first, not in the lowest, Landau level. This does not occur 
for the cases shown here, but a simple estimate shows that it occurs if the anomalous moment and the charge have 
opposite sign and |k| becomes of the order of or larger than the modulus of the normal magnetic moment ^j^. 

One can check that for large magnetic fields the onset lines are well approximated by the curves /i» = m» — kB, 
which turn out to become straight lines for large B. However, we need to remember that we have used an effective 
approach for the anomalous magnetic moment. As a consequence, we should not trust our results for kB becoming of 
the order of or larger than m*. In that case, as we see for instance from eq. (|61[) . M l/= o iS= +i would become negative 
(for k > 0). In the case of proton matter we find kB ~ 0.4m» at B = 5 x 10 19 G, which suggests that our approach 
is still reliable in the plotted range. We expect that in a full treatment the onset lines saturate for asymptotic values 
of B. This is the case for k = 0, where our approach can be used for arbitrarily large B. 



C. Neutral baryons 



For the case of neutral baryons we may think of pure neutron matter with q = and k = — 1.91 fj,N- In this case 
there is no Landau quantization and the single-baryon excitations become 



£k. s = \ k\ 



k 2 , — skB 



(65) 



where kj_ is the momentum perpendicular to the magnetic field. Replacing the momentum integral in the nucleonic 
pressure (|59)l . 



d 3 k 



—y 

2tt 2 ^ 



s=± 



fOO 

dkiki I dk\\ 



(66) 
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Hb [MeV] B [10 iy G] 

FIG. 4: Same as fig. [3] but for neutral baryons, q — 0. Due to the symmetry under k — ¥ —K, we can restrict ourselves to 
negative k. The density along the onset in the right panel is shown for pure neutron matter. We have plotted the total baryon 
density (solid) and the contributions from the s = +1 and s = —1 states (dashed). We see that at B > 0.4 x 10 19 G the system 
is fully polarized in the s = — 1 state. 



and performing the integrals at zero temperature yields 

6(m* + SK B — m, 



Ft, 



E 

s=± 



+ 



2tt 2 
skB^I 



k 2 F a Af 2 skBM s 



12 8 



M * , j r3 ( M s , skB \ i M* + k F. t 

arccos h M, In 

6 /i* s \ 8 6 J M a 



(67) 



with M s = m* — sk-B and the longitudinal Fermi momentum kp s = y^u| — M 2 . The baryon and scalar densities 
become 
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We observe that the system is symmetric under k — > — k because after a sign change of k the contributions from 
s = +1 and s = — 1 have simply exchanged their roles. Analogously to the charged case we can now compute the 
baryon onset. The results are shown in fig. [4] For large magnetic fields, the onset curves can be approximated by 
the simple straight lines /ib^tob — \ k\B. In this regime the system is fully polarized, i.e., the baryons are all in the 
s = — 1 (s = +1) state if k < (n > 0), as can be seen in the right panel of the figure. Having in mind the symmetry 
under k — > — K it is instructive to go back to the results for charged baryons. In the left panel of fig. [21 the onset lines 
for k = +(1n and k = — /ijy are obviously far from identical. This is a LLL effect: in the LLL the baryons with q > 
become heavier for k < and lighter for k > 0. If we let q — > the higher Landau levels become more and more 
important, and the two curves indeed approach each other and become identical for q = 0, see onset line for n = — /ijv 
in fig. |U Now the baryons in the lowest energy state become lighter for either sign of n. 



D. Binding energy 

Let us now compare the baryon onset in the holographic model, fig. [2 with the ones in the Walecka model for 
charged and neutral baryons, figs. |3] and |4l The first simple observation is that in the holographic case the larger the 
magnetic field, the larger the energy needed to create baryons. This seems to be in contrast to the results from the 
previous two subsections which suggest that a magnetic field tends to make nuclear matter energetically less costly. 
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FIG. 5: Solid lines: binding energy Eq along the baryon onset for charged baryons with charge q = +e and various anomalous 
magnetic moments k, corresponding to the four cases shown in fig. [3] Dashed lines: binding energies for baryons with the 
same values for k, but q = 0. They approximate the curves for the charged baryons at small B where many Landau levels are 
occupied. The arrow at Eq = —56.3 MeV indicates the asymptotic value of the binding energy for k = 0, see eq. (|70[) (for k/0 
our approach does not allow for arbitrarily large magnetic fields). 



In order to compare the holographic with the Walecka result in a sensible way, we need to isolate the effect of the 
binding energy Eq- 

In the absence of a magnetic field the binding energy is given by Eq = ^—niB, i-e., by the energy per baryon = 
relative to its mass ms- To be precise, by mass we mean the energy that is needed to put a single, non-interacting 
baryon into the lowest single-particle state of the system. To generalize this concept to finite magnetic fields, let us first 
consider charged baryons. In this case, the single-particle ground state energy without nucleon-nucleon interaction is 
tub — nB&gnq according to eq. (|61j) . Thus we define 

E = — (tub - kB sgn q) , (69) 

n B 

where, at the onset, ^- = \Ib- We plot Eq as a function of the magnetic field in fig. [5] The plot shows that in the 
regime where higher Landau levels are occupied, the binding energy depends strongly on the anomalous magnetic 
moment, while it varies very little with k in the LLL regime. For k = we can take the limit of asymptotically large 
magnetic fields and find that the onset approaches the line /x» = to*, hence with the definitions of /x* and m* below 
eq. (|59|) we have 

Eq(k = 0,B -> oo) = g^ujQ - g a o ~ -56.3 MeV. (70) 

Here, ujq and a arc complicated functions of the parameters of the model. Nevertheless, this expression is very 
instructive since it shows the effect of the meson condensates in a very transparent way: E is negative because the 
scalar meson condensate, responsible for the attractive interaction, becomes sufficiently large compared to the vector 
meson condensate, which is responsible for the repulsive interaction. There is no such simple expression for B = 0. 
Interestingly, Eq(B — > oo) is nonvanishing only in the presence of scalar self-interactions. For asymptotically large 
B, one can show analytically that after setting b = c = in eq. (|58[) and (|60ap we would obtain \ib = tib and thus 
So =0. 

As a result of this discussion we conclude that the magnetic field has two effects: it changes the mass and the 
binding energy For charged baryons whose charge and anomalous magnetic moment have the same sign, both effects 
work in favor of creating baryonic matter, at least for sufficiently large magnetic fields: first, they decrease the mass 
itib — nBsgnq, and second, they lead to a larger binding energy \Eq\. At smaller magnetic field and/or different signs 
of q and k, things are more complicated and can be read off from figs. [3] and [5j For neutral baryons, a magnetic 
field always decreases the mass, niB — \k\B, and baryonic matter is always favored by a nonzero magnetic field. If we 
compare charged with neutral baryons at the same anomalous magnetic moment, for instance the curves for k = +/-J>n 
in figs. [3] and |4l we see that neutral baryonic matter is favored a bit less than charged one. The reason is that for 
neutral baryons jii'ol becomes smaller for large magnetic fields, as fig. [S] demonstrates. 
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E. Large N c 

How does this picture change when we go to the large-7V c limit? One way to answer this question is to si mply 
rescale the parameters of the model with the appropriate powers of N c , as suggested from large- arguments [6711681]. 
This has been done without magnetic field in ref. |53|. While the rescaling of most parameters is unambiguous, one 
has to make a decision about the generalization of the scalar meson to large values of N c . In the traditional quark- 
antiquark picture, its mass would scale like N®. However, we may take the heaviness of the lightest scalar meson in 
the Sakai-Sugimoto model as a hint that this picture is incorrect, see also refs. (6^-tZll- The alternative that we shall 
consider here is that the lightest scalar meson is a tetraquark state [72| . Let us denote the generalization to arbitrary 
values of N c of this tetraquark state by v. This state is composed of N c — 1 quarks and N c — 1 antiquarks such that its 
mass scales like m x ~ 2(N C — 1) ~ N c |73|. There are other possibilities for the nature of the scalar meson which we 
do not consider here. The results of ref. [53|, however, suggest that the main conclusions of our followin g d iscussion 
do not depend on the detailed nature of this controversial meson state. Consequently, we rescale [53l. l67l. l68j 

m B ,m x ,K,q ~ N c , m u ,g x ~ N° , g u ~ N? 2 . (71) 

(The iVc-dependence of K has been computed in ref. [74j : note that we also assume the charge q to scale with N c ). 
The 7V c -dependence of the self-interactions of x is given by b ~ e~ Nc , c ~ N c ; this is suggested by the arguments 
explained in ref. [53j . 

Let us first briefly discuss the scenario without magnetic field. In this case the stationarity equations of the Walecka 
model (|60p imply that the Fermi momentum kp = y//f — mf scales like N®, although both /i* and m* arc proportional 
to N c . Therefore, in the large- N c limit, one can expand the pressure for small kp <C tti*, /x» (which, as an aside, yields 
P ~ N c ). The resulting equation P = for the baryon onset then becomes very simple and yields, together with 
the stationarity equations the trivial solution m* = ra^, /i* = fig and E Q = 0, i.e., the baryon onset becomes a 
second-order transition. This result is in fact obtained for all N c > 4 [53j . 

Now we switch on the magnetic field. Our numerical results show that, not surprisingly, the binding energy remains 
zero, and as a consequence the baryon onset curves are simply straight lines in the P>-\ib plane for all N c > 4: for 
charged baryons they are given by B = mB ~^ B ; while for neutral baryons B = mB \~\ lB ■ Note the slight subtlety 
related to the order of limits q — >• and N c — s- oo: consider the onset lines of charged baryons at N c = 3 for k = ±/Uat, 
see fig. [3] These are two very different lines. As discussed above, these lines merge for q — ¥ 0. They remain on top 
of each other if now we let N c —> oo. On the other hand, if we first take the limit N c — s- oo they remain separated 
and become linear with opposite slopes. Now letting q — > does not change the result. Hence the symmetry of the 
onset magnetic field under n — > — k if we take the large- N c limit of our neutral baryons and the aritz-symmetry for 
the charged baryons. 

With these preparations, what can we learn from our holographic result? We know that there is no binding energy 
in this case, and the onset line simply indicates the B-dependent baryon mass (at least without meson supercurrcnt, 
otherwise the situation is more complicated, as discussed in sec. HIIBj) . This mass gets larger with increasing magnetic 
field. In the Walecka model, the only case with increasing mass is the case of charged baryons with q and k having 
opposite sign. Since in the present context a nonzero electric charge is equivalent to the existence of Landau levels, 
it is interesting to ask whether there is any Landau level structure for our holographic baryons. We know that in the 
chirally restored phase of the Sakai-Sugimoto model there is indeed a phase transition which has similar properties as 
a transition into the LLL [27], E^J . In terms of the solution to the equations of motion, the "LLL" phase corresponds to 
a trivial solution, = oo in the notation of ref. [27j |. This solution always exists, but in some regions of the parameter 
space is disfavored compared to a nontrivial solution, interpreted as a phase of higher Landau levels. In the baryonic 
phase, there is also a trivial solution, namely j/oo = 0, see eqs. (|35[) . However, this solution is unphysical because 
it leads to an infinite baryon density and thus infinite free energy Although we find several nontrivial solutions in 
certain parameter regions, the solution which is continuously connected to the solution at the onset is always the 
energetically preferred one. Hence there is no phase transition within the baryonic phase. This is consistent with the 
apparent bosonic nature of the holographic large- N c baryons (see sec. IIII A|) because for a Bose condensate at zero 
temperature we do not expect de Haas-van Alphen oscillations. 

Besides the overall tendency of the B-dependent mass, how about its linear behavior seen in the large- N c limit of 
the Walecka model? The holographic results show an approximately linear behavior only for intermediate magnetic 
fields. For large fields, a comparison to the Walecka results would only be sensible if the latter included a more 
elaborate treatment of the anomalous magnetic moment. But also at small magnetic fields our holographic onset line 
differs from the linear behavior, in fact we find /x on set = + const x b 2 + . . .. It would be interesting to compute 
the onset in other field-theoretical approaches. After all, the Walecka model does not know about dynamical chiral 
symmetry breaking and thus cannot be expected to show effects from MC and the meson supercurrent, which, as 
mentioned above, seem to be the driving forces for the independence of our effective holographic baryon mass. On 
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the other hand, also our holographic approach should be refined for a more meaningful comparison. In particular, a 
generalization to two flavors is necessary to describe realistic nuclear matter. 



V. PHASE DIAGRAM WITH HOLOGRAPHIC BARYONIC MATTER 



We now come back to the Sakai-Sugimoto model. In this section we are mainly interested in the chiral phase 
transition in the presence of baryonic matter. This requires the inclusion of chirally symmetric matter into our 
analysis. Since we consider both broken and symmetric phases in the deconfined background, the metric does not 
change. However, the embedding of the flavor branes into the background is different. In the Sakai-Sugimoto model, 
chirally symmetric matter is described by a configuration where the D8- and D8-brancs are disconnected and have 
a trivial embedding in the background geometry, x'^(u) = 0. As in the chirally broken phase, we shall work in 
the approximation f(u) ~ 1. In contrast to the broken phase (see discussion at the beginning of sec. Ill B[) this 
approximation corresponds unambiguously to the zero-temperature limit. The reason is that, since the flavor branes 

3 

do not connect, the holographic direction on these branes is parameterized by u e [ut, oo] ■ Hence, f(u) = 1 — -J can 
only be approximated by 1 if ut —> 0, i.e., T — > (see sec. 4 of ref. (27j as well as ref. [4|| for the full temperature- 
dependent treatment and the limit T — > 0). 

Our main result of this section will thus be the phase diagram in the b-\x plane at zero temperature. For our purpose 
this is not a severe restriction for the following reasons. Judging from our previous results without baryonic matter, 
the most interesting physics is observed at T = 0. For instance, IMC is most pronounced at T = 0. The main effect 
of nonzero temperature is to disfavor the chirally broken phase, and the non-monotonic behavior of the T = chiral 
phase transition line is weakened, but not altered in a nontrivial way. Moreover, since our main results concern matter 
at large chemical potential, the only application will be in neutron stars where the temperature is negligibly small 
compared to the chemical potential. 



A. Results 



For the phase diagram, we have to compare the free energies of three phases, corresponding to three different 
embeddings of the flavor branes, V, U, and ||, 

fi v f°° 3/2 u 5 + b 2 u 2 - \[(bu c + d) 2 -n 2 } | n 4 u c fin 

N L s/giu) 6 2 ' (72a) 

Ou = r ,,„„,, u^b 2 u 2 -^ { ul + b 2 u 2 ) , n (?2b) 

N J™ ^Ju% + b 2 U$ - (ug + b 2 ul) - (u 3 - U 3)_jL.( u 5 + 6 2 M 2) 2 

O r°° „,5 I 1.2 2 | 1 (3&/-Q 2 

iuu^ +bU +2 ^, ^ (72c) 



U Jo / u S + b 2 u 5 +u 3j3b^ 

y sinh^ z c 



The first free energy is taken from sec. Ill Cl and describes baryonic matter, i.e., the flavor branes are connected and 
have a cusp in their profile at the location of the pointlikc D4-brancs u = u c . The second describes mesonic matter 
which we have worked out in ref. [27j |, In this case, the flavor branes are also connected but have a smooth profile 
everywhere, in particular at tip of the branes u = uq. The quantities Uq and rj, appearing in the free energy, have 
to be determined from two coupled algebraic equations. At the baryon onset, we have u c = uq and rii = 0. Using 
the explicit form of the equations for uq and rj [27|, this yields {bu c + d) 2 — n\ = (u^ + 6 2 Uq), which shows that 



at the onset f2 v = ^u- The third free energy is also taken from ref. |27J and describes chirally restored matter, i.e., 
the case of disconnected flavor branes. Here we have one algebraic equation which has to be solved numerically for 
the quantity Zoo . Since different solutions constitute the state of lowest free energy in different regions of the phase 
diagram, there is a first-order phase transition within the symmetric phase which can be associated with a transition 
into the LLL. While the phase of higher Landau levels ("hLL") corresponds to a nontrivial solution, the "LLL" phase 
has the trivial solution z,^ = oo, which leads to a simple expression for the free energy, 

= r du Vu 5 + b 2 u 2 - ^ . (73) 
■A/ Jo 2 
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FIG. 6: Zero-temperature phase diagram in the plane of magnetic field b and quark chemical potential /i without (left) and with 
(right) baryonic matter in the deconfined geometry of the Sakai-Sugimoto model. In shaded areas the ground state is chirally 
symmetric, while in unshaded areas chiral symmetry is spontaneously broken. Solid lines are first-order phase transitions, the 
dashed line is the baryon onset. The two different chirally symmetric phases are reminiscent of phases where only the lowest 
Landau level ("LLL") and where higher Landau levels ("hLL") are occupied. Recall that b and [i are dimensionless quantities 
that involve dimensionful parameters of the model. As a rough estimate, b — 0.1 corresponds to \q\B ~ 5 x 10 18 G. (This 
estimate is obtained by using the simple fit of ref. [53] which involves the assumption that the 6 = chiral phase transition 
without baryonic matter at fj, « 0.44 occurs at a quark chemical potential of 400 MeV.) 
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FIG. 7: Left panel: location of the tip of the connected branes u c (denoted uo in the mesonic phase) as a function of fj, for 
three values of the magnetic field, i.e., along three horizontal cuts through the phase diagram in the right panel of fig. [6] 
One representative for each of the three qualitatively different cases is shown: baryon onset (small magnetic fields); baryon 
onset followed by chiral phase transition (intermediate magnetic fields); chiral phase transition (large magnetic fields). As fig. 
[6] shows, for chemical potentials beyond the scale of the plot the system reenters the chirally broken phase, resulting in an 
additional first-order phase transition in all three cases. Right panel: baryon number density n along the same cuts through 
the phase diagram. Through the axial anomaly the meson supercurrent produces a nonzero n also in the mesonic phase for 
b, n > 0. 



Obviously, this free energy, like all three expressions in eqs. (|72|). is divergent. After subtracting the vacuum free 
energy flu(p = 0) all expressions become finite. Equivalently, we can simply take the pairwise difference of free 
energies to find the ground state at any given b and \i. The resulting phase diagram of this numerical calculation is 
shown in fig. [5] 

For comparison, we have also included the phase diagram without baryonic matter from ref. {27} . In fig. [7] we 
present the location u c of the tip of the connected branes and the baryon density along lines of constant magnetic 
field. The interpretation of the former deserves a comment. In the mesonic phase, u c is usually identified with the 
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constituent quark mass since it is the distance between the color and flavor brancs. This is still true in the baryonic 
phase, but there it has a second meaning because Nc ™ c is the baryon mass. (Note that the factor -| originates from 
the geometry of the model and has nothing to do with the number of colors.) 



B. Observations 



The main observations can be summarized as follows. 

• The effect of baryonic matter at small magnetic fields is dramatic: it prevents the system from restoring chiral 
symmetry for any chemical potential. In doing so, it completely expels the "hLL" phase from the phase diagram, 
such that the only surviving chirally restored phase is the "LLL" phase. 

• The baryon onset line terminates at a critical endpoint (n,b) ~ (0.26,0.25) where it intersects the chiral phase 
transition line. As a consequence, at sufficiently large magnetic hclds the mesonic phase is superseded by quark 
matter before baryonic matter can even be created. 

• In the presence of baryonic matter, IMC plays an even more prominent role in the phase diagram: at any given 
chemical potential /i > 0.26 a sufficiently large magnetic field induces chiral symmetry restoration. 

Let us briefly comment on the first of these observations whose 6 = version was already made in ref. [47j . The fact 
that chiral symmetry remains broken along the entire \i axis is an apparently puzzling result: in view of the similarity 
of the present non-antipodal version of the Sakai-Sugimoto model and the NJL model, one might have expected chiral 
symmetry to be intact for sufficiently large chemical potentials. This was true before including baryonic matter, 
in which case the entire phase structure with magnetic field was in amazing agreement with NJL calculations f27j . 
However, in the presence of baryonic matter it is questionable to expect this agreement to persist for the following 
reason. In related NJL studies PjJJUl Ell baryon number is a rather simple concept: the only degrees of freedom are 
quarks whose masses acquire a contribution from the chiral condensate. When the chemical potential is larger than 
this constituent quark mass, a nonzero quark number is generated. Then, baryon number is simply this quark number 
divided by N c . This is not true in our present approach. Holographic baryons are different from a mere collection of 
N c quarks. This is clear from their construction and can easily be seen from the difference between the constituent 
quark mass in the mesonic phase u c and the baryon mass |" c . Whether we should therefore call our quarks confined 
is debatable, but it is certainly a qualitative difference to the NJL model. 

To interpret the lack of chiral symmetry restoration, the following intriguing property of our baryonic matter at 
asymptotically large chemical potentials may be helpful. For simplicity, let us restrict ourselves to b = 0. In this case, 
as demonstrated in appendix [Cl the free energy of the baryonic phase approaches the one of the chirally restored 
phase, 

where p = T (•jjj) T (|) /^/n (see fig. |8]in appendix lO for a plot of the full numerical result). Therefore, at asymptot- 
ically large p, baryonic matter and quark matter become thermodynamically indistinguishable. One may speculate 
that this is a consequence of the pointlike nature of our baryons: due to this property, our baryons can only overlap 
at infinite density. This suggests that the expected transition to quark matter at finite fi is shifted to (J, = oc (which, 
curiously, is undone by a sufficiently large magnetic field). It would be interesting to see whether this is different for 
Nf > 1, where baryons can be described by non-singular instantons. 



VI. SUMMARY AND OUTLOOK 



We have discussed baryonic matter in a magnetic field, using the deconfined geometry of the Sakai-Sugimoto model 
where baryons are introduced by D4-brancs wrapped on the internal S* 4 . Our main focus has been the onset of 
baryonic matter and its effect on the chiral phase transition at zero temperature and finite chemical potential. 

The critical chemical potential for the onset of holographic baryonic matter increases monotonically with the 
magnetic field, saturating at a finite value for asymptotically large magnetic fields. For subcritical chemical potentials 
the system is in the mesonic phase with a meson supercurrent in the direction of the magnetic field. Because of the 
axial anomaly, the system has nonvanishing baryon number for all nonzero magnetic fields and chemical potentials. 
Due to the presence of the supercurrent the baryon onset is significantly "delayed" to larger chemical potentials. In 
contrast to real-world baryonic matter, the onset is a second-order phase transition. Besides the holographic study, 
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we have computed the onset also in the Walecka model, a hcld-thcorctical model for dense nuclear matter where the 
onset is a first-order transition. Within this model, we have demonstrated that the onset depends strongly on the 
electric charge and the anomalous magnetic moment of the baryons. In most cases, however, the magnetic field favors 
baryonic matter because it decreases the baryon mass, although for neutral baryons it also decreases the (modulus of 
the) binding energy. In the holographic calculation there is no binding energy and the baryon mass is increased by 
the magnetic field. Our results indicate that this increase is closely related to the effect of MC. 

Without baryonic matter, chiral symmetry is restored for any given magnetic field at some sufficiently large chemical 
potential (27j . With baryonic matter, this is only true for sufficiently large magnetic fields, where we have found that 
baryons play no role. For small magnetic fields there is a transition from mesonic to baryonic matter but no subsequent 
transition to quark matter. This enforces the unusual effect of IMC: although the magnetic field typically enhances 
a chiral condensate, it now restores chiral symmetry for any given chemical potential larger than, roughly speaking, 
the one where the baryon onset line and the chiral phase transition line meet in a critical endpoint sec fig. [B] 

As our comparison with the Walecka model has shown, one has to be very careful with drawing conclusions from 
the holographic results for the QCD phase structure and the interior of neutron stars, where the extreme conditions 
studied here might be realized. Most importantly, our holographic calculation has been restricted to the N c — > oo 
limit, and there are indications that large-iV c nuclear matter is very different from real nuclear matter. Therefore, 
generalizations to finite N c would be very interesting, as it has been done for instance in related D3/D7 models [75| . 
On the other hand, in certain aspects our approach seems to be more realistic than widely used models of dense 
matter. In particular, it is interesting to compare our results to the NJL model. In our previous work [27J we have 
pointed out an amazing agreement with corresponding NJL phase diagrams [IJ G3- This agreement is lost after 
including baryonic matter. This is understandable since in the NJL model there arc no baryons; baryon number is 
only generated by deconfined quarks. In the Sakai-Sugimoto model, however, baryons are clearly distinct from a set 
of N c deconfined quarks. 

It would be interesting to extend our study to more than one flavor in order to describe more realistic baryonic 
matter. Moreover, one might include superfluidity of nuclear matter, as suggested for the present model in ref. [76j . 
In the Sakai-Sugimoto model as well as in AdS/QCD approaches, it has also been suggested that - even for vanishing 
ma gne tic field - the ground state breaks rotational and/or translational invariance at sufficiently large baryon densities 
[77H8G ] . It would be interesting to study the effect of such phases on our phase diagram. We should also keep in mind 
that the quark matter phase considered here is not in a state expected from QCD. As indicated in sec. II A[ cold and 
sufficiently dense quark matter is a color superconductor in the CFL phase. This phase also breaks chiral symmetry 
at asymptotically large fi, like our holographic baryons, but the mechanism is very different and heavily relies on the 
fact that N c = 3. 
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Appendix A: Solving the equations of motion 

In this appendix we solve the equations of motion ([TO]) . To some extent we can use the same method as in ref. [27j 
for the case without baryons. It is useful to rewrite the equations as 

(Ala) 
(Alb) 
(Ale) 



We integrate eq. (|Ala[) to obtain 
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where the integration constant k can be determined from evaluating this equation at u = u c , 



d±^£], (A3) 



u 

"3 V~ 



where a3(u c ) = and ao(u c ) = ±^ has been used. Using eqs. (|Ala[) , (jAlbp . (|A2[) , and (|A3|) . we compute 
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Inserting this expression into the original equations of motion for ao and 03 (obtained from taking the derivative with 
respect to u on both sides of eqs. (|19a[) and (|19bj) ). one shows that these two equations are equivalent to 



d y a = a 3 (j/) + ^i, (A5a) 
d 



d y a 3 = a (y) + — , (A5b) 
where the new variable y fulfils 

y' = -™ , (A6) 

V9W 

with g(u) defined in eq. (|29|) . Now, inserting a' 3 = y'd y a 3 = y' (ao + ^) into eq. (jAlcft yields the solution for the 
embedding function, 

f u dv 

xAu) = k — - — . (A7) 

Ju c V*/ 2 ^V) 

The system of coupled equations (|A5|) for a and a 3 can easily be solved. With y(u c ) = we have the boundary 
conditions ao(y = 0) = ±3^, a 3 (y = 0) = 0, and the solutions become 

ao(y) = — [(d ± 6u c ) coshy + n 4 sinhy - d] , (A8a) 

a 3(y) = — [(d±bu c ) sink y + n4(coshy- 1)] . (A8b) 
36 

With j/oc = y(u = 00) the boundary conditions at the holographic boundary ao(yoo) = /1, 03(2/00) = J can be used to 
determine d and n 4 , 

d = ~2 M ~ 2 (, T ~ J rinhy, J ' (A9a) 

3, . 3 / ^ it c \ 1 + cosh?/ 00 
" 4 = -2 &J+ 2H /i=F T) sinh^ ■ (A9b) 

These expressions are valid for any j. From the main text we know that the condition that j minimize the free energy 
is equivalent to d = — Inserting this value into eq. (|A9a|) yields the supercurrent given in the main text, eq. (|36p . 
and, inserting the result for j into eq. (|A9b[) . the density 714 given in eq. (|34[) . The results for d and 714 can then be 
inserted into eqs. (|A8|) to obtain eqs. (f3Tj) . 

Comparing eq. ([22]) with eq. (|A9b[) we obtain the total density 

n=-b(n=f Uc ) 1 + cosh y°° (A10) 
2 V 3 / sinhj/oo 

With the help of eqs. (|A9[) and fc from eq. (jBDj) the solutions for X4, ao, 0,3 in eqs. (|A7|) and (|A8[) can be written in 
terms of u c and j/oc (and and 6, which are externally fixed variables). The remaining two quantities, u c and yoo, 
need to be determined numerically from the following two coupled equations: firstly, the defining equation for y^, 

f°° u 3 / 2 du 

Voo = Zb !L=™, (All) 

and secondly, the equation that fixes the asymptotic separation of the flavor branes, Ijl = x 4 (oo) — a; 4 (w c ), which 
reads 

£ , [°° du IK oX 

- = k ■ . A12) 

2 J Uc U^^gju) 
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Appendix B: Approximations close to the baryon onset 
1. Vanishing B 

In this appendix we derive eq. (|47|) . i.e., the behavior of the baryon density close to the onset at vanishing 
magnetic field. For b = and with = Sby 1 ^, cqs. (|35p can be written as 



L 3 ' 2 du 



v^-i + R^ 3 -!)' 



(Ill 



9 ^ 5J 1 u^^- 1 + ^(^-1) 



(Bla) 



(Bib) 



For chemical potentials close to, but above, the onset chemical potential ^Q Ilsct 
the ansatz 



Voo — Voo, onset "1" ' 



with a, (3 to be determined, and 



Inserting this ansatz into eq. (|44[) yields 



fi-m g . 



m® the numerical results suggest 



(B2) 
(B3) 



a 



3?;° 

^iJoc. onset 



(B4) 



To compute a it is sufficient to consider eq. (|Blb[) . where j3 does not appear. We are interested in the terms linear 
in e. The linear term on the left-hand side is easily obtained. On the right-hand side, we subtract the constant term 
and neglect the quadratic term in the square root in front of the integral. If our ansatz is correct, the integral then 
must yield a linear term but, as written, also yields terms of higher order in e, 



at 



4«o„ S ct) 1/2 



du 



j,3/2yJ U 8 -l+V 2 e 2 ( U 3 _ |) 



Here we have abbreviated 



Q„0 ( v \5/2 

oo. onset V c, onset ) 



(B5) 



(B6) 



Let us for the following arguments denote the integral in eq. (|B5j) by X. One can check numerically that I is indeed 
linear in e for small e. Obviously, we cannot proceed by naively expanding the integrand in e since this procedure 
would miss the linear term. Instead, we employ the following trick. Neglecting higher order terms, I divided by 
e should not depend on e anymore, i.e., J^j ~ 0. We evaluate this equation by first rewriting the derivative with 
respect to e as a derivative with respect to u. 



d 



1 



2ev 2 (u 3 



d 



1 



de y/u*-l + « 2 e 2 (u 3 - |) " 2 ( 8u5 + 3y2e2 ) du y/u*-l + «V (u 3 - I) 
Then, with partial integration we obtain 



(B7) 



I 
r 



"12 



du 



1 



d 



^ U 8_ 1+U 2 C 2 (U 3_|) du uV^ + ?,V^) 



(B8) 



where the first term on the right-hand side is the boundary term. Now the left-hand side and the boundary term are 
constant in e while the second term on the right-hand side is linear in e and can thus be dropped (the integral is finite 
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and yields a term constant in e, but no e 1 term). Consequently, we have arrived at the very simple result I : 
which can be confirmed numerically and which we insert into eq. (|B5[) . The result is 



q?;° (n° VP - 1 

^tfoo.onsctV^c, onset/ *- ± 



(B9) 



With the same trick eq. (|Bla[) can be evaluated to obtain /3. Here we are only interested in the baryon density for 
which we insert eq. (|B9[) into eq. (|B4[) . With it° onset and y^ onset from eq. (|45|) we obtain the final result given in eq. 
(|4"T)) in the main text. 



2. Asymptotically large B 



Here we derive eq. (|52p . i.e. the behavior of the usual baryon density 714 close to the onset at asymptotically large 
magnetic fields. In this case, due to the supercurrent, the onset does not occur at the baryon mass but at twice the 



baryon mass, 



2m™, where 



t /3 denotes the (dimensionless) constituent quark mass in a baryon 



at b — > 00. As in the previous appendix, we first need to compute j/oo and u c . With y^ — > 00 for 6 
from eq. (J31J that 
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Y 1 



2iC 



we conclude 



(BIO) 



and thus we can write the second of eqs. (j3"5")) as 



1/2 



1 - 



36w= 



da 



.3/2 



(Bll) 



1 



This single equation can now be used to determine the behavior of at the onset. Again, as for b = 0, the numerical 
result suggests the ansatz 



ae , 



e = fx- 2m" 



with 5 to be determined. Inserting eq. (|B12|) into eq. (|B10|) yields 

5 = 1- 



36 2a 
, 4 c-ve, v = l--. 

With the help of these expressions for n|° and we see that to zeroth order in e eq. (|Blip implies 

leTrfrd)' 



16tt 

isct p2 



r(A) 



(B12) 



(B13) 



(B14) 



from which the onset chemical potential in eq. (|5 1 1) is obtained. It is obvious that the left-hand side of eq. (|B11[) has 
a linear term in e. For the right-hand side, however, we need to apply a similar trick as in the previous appendix to 
extract the linear term. We write the linear terms of eq. (jBllj) as 



4(u= 



\2 
b 1 



0_ 

8e 



du 



u 5 - 1 



4R 



(B15) 



e=0 



Now, analogously to eq. (|B7|) we rewrite the differentiation with respect to e by a differentiation with respect to u 
and compute the integral via partial integration. Then, for e = only the boundary term survives, and we can easily 
compute the result for a, 



12 



8 - 15« onS et) 1/2 ^ 

Inserting this into eq. (|B12[) and the result into eq. (|B10[) yields n|° as given in eq. ([52")) in the main text. 



(B16) 
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b=0 



0.5 - 



0.0 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
20 40 60 80 100 

FIG. 8: Ratio of the pressures of baryonic matter P v and quark matter P|| at vanishing magnetic field as a function of fi. The 
plot shows that baryonic matter always has larger pressure, i.e., chiral symmetry is never restored, but at asymptotically large 
fj, baryonic and quark matter become indistinguishable. (Note the huge /i scale compared to the phase diagrams in fig. [6]) 



Appendix C: Asymptotics at large chemical potential 

In this appendix we compute the behavior of the baryonic phase at asymptotically large fi and b = 0. The numerical 
evaluation shows that j/oo oc /i~ 3 ^ 2 , u c oc fi° for fj, — > oo. Hence we make the ansatz 

y°° - -^p- » ( C1 ) 

where p is a number which we shall now determine. We have extracted the linear behavior in b which is necessary to 
take the 6 — s- limit. Inserting this ansatz into eq. (|34|) yields 

5/2 

(1 > 



n 4 ~ iq , (C2) 
and thus the equation for ([55]) becomes 



9 /.I- 



with the new integration variable x = . Now we can approximate the lower boundary by and neglect the term 
oc /j~ 3 in the denominator of the integrand. Then, performing the resulting integral, we see that u c drops out of the 
equation and we obtain 

r (±) r (£) 

= \wl± ^ . (C4) 



Next we can compute the free energy. From the general expression (|40| we find with the above asymptotic relations, 
the same change of integration variable, and recalling that n — at b = 0, 



Here we have performed the integral explicitly with a cutoff A at the upper boundary and noticed that the asymptotic 
b = result of the baryonic phase is exactly the same as the full = result for the "hLL" phase with restored chiral 
symmetry, see appendix D of ref. [27j • At non-asymptotic values of (U, the free energies of the two phases differ. In 
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fig. [8] we plot the ratio of the two corresponding pressures P = — Q, obtained from the numerical result. 
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